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ABSTRACT 

The BABAR experiment has discovered a new narrow state, -D*j(2317) + , near 2.32 GeV/c 2 
in the inclusive D+ir invariant mass distribution from e + e~ annihilation data at 
energies near 10.6 GeV |T]. The same experiment has also shown evidence for struc- 
ture in the 2.46 GeV/c 2 region in the D*(2112) + 7r° mass spectrum. These discover- 
ies have triggered several experiments in a search for new states coupled to the -D+ 
meson which confirmed the existence of _D*j(2317) + together with D s j(2A58) + — > 
D*(2112) + 7T° both in inclusive e + e~ annihilation and in B decays. These two new 
states are difficult to explain in terms of potential models. 
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1 Introduction 



Experimental information on the spectrum of the cs meson states is limited. The 
ground state, the meson, is well-established, as is the 3 Si ground state, the 
D*(2112)+. Only two other cs states are listed in the last PDG edition [2j. The 
P s i(2536) + has been detected in its D*K decay mode and analysis of the D* decay 
angular distribution prefers J p = 1 + . The _D*j (2573) + was discovered in its D°K + 
decay mode and so has natural spin-parity. The assignment J p = 2 + is consistent 
with the data, but is not established. 

The spectroscopy of cs states is simple in the limit of large charm-quark 
mass [3 Ej. In that limit, the total angular momentum j = I + s of the light 
quark, obtained by summing its orbital and spin angular momenta, is conserved. 
The P-wave states, all of which have positive parity, then have j = 3/2 or j = 
1/2. Combined with the spin of the heavy quark, the former gives total angular 
momentum J = 2 and J = 1, while the latter gives J = 1 and J = 0. The J p = 2 + 
and J p = 1 + members of the j = 3/2 doublet are expected to have small width [Sj, 
and are identified with the D*j(2573) + and -D s i(2536) + , respectively, although the 
latter may include a small admixture of the j = 1/2, J p = 1 + state. Theoretical 
models typically predict masses between 2.4 and 2.6 GeV/c 2 for the remaining two 
states 13 Ej, both of which should decay by kaon emission. They would be 
expected to have large widths j3 E] and hence should be difficult to detect. 

The experimental and theoretical status of the P-wave cs states thus can 
be summarized by stating that experiment has provided good candidates for the 
two states that theory predicts should be readily observable, but has no candidates 
for the two states that should be difficult to observe because of their large predicted 
widths. 

2 The BABAR Experiment. 

The BABAR detector (at the PEP-II asymmetric-energy e + e~ storage ring with 
center-of-mass energy near 10.6 GeV) is a general purpose, solenoidal, magnetic 
spectrometer and is described in detail elsewhere [SJ. The data sample use in this 
analysis corresponds to an integrated luminosity of 91 fb~ . 
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3 Dfir events selection. 

The objective of this analysis is to investigate the inclusively-produced Dfn mass 
spectrum by combining charged particles corresponding to the decay — ► K + K~tt + 
with 7r° candidates reconstructed from a pair of photons and performing a one- 
constraint fit to the 7T° mass. A given event may yield several acceptable ir° candi- 
dates retaining only those candidates for which neither photon belongs to another 
acceptable tt° candidate. To reduce combinatorial background from the continuum 
and eliminate background from 5-meson decay, each K + K~tt + 7t° candidate was 
required to have a center-of-mass momentum p* greater than 2.5 GeV/c. 

The upper histogram in Fig.^a) shows the K + K~ir + mass distribution for 
all candidates. Clear peaks corresponding to D + and Df mesons are seen. To reduce 
the background further, only those candidates with K + K~ mass within 10 MeV/c 2 
of the 0(1020) mass or with K~ir + mass within 50 MeV/c 2 of the if* (892) mass are 
retained. The decay products of the vector particles 0(1020) and if* (892) exhibit 
the expected cos 2 6^ behavior required by conservation of angular momentum, where 
6h is the helicity angle. The signal-to-background ratio is further improved by 
requiring |cos#/J > 0.5. The lower histogram of Fig. ^a) shows the net effect of 
these additional selection criteria. The Df signal and sideband regions are shaded. 
The signal peak, consisting of approximately 80,000 events, is centered at a mass 
of (1967.20 ± 0.03) MeV/c 2 (statistical error only). 

Figure dh) shows the mass distribution for all two-photon combinations 
associated with the selected events. The 7r° signal and sideband regions are shaded. 
Candidates in the signal region of Fig. ^a) are combined with the mass- 
constrained 7T° candidates to yield the mass distribution of Fig. ^c). A clear, nar- 
row signal at a mass near 2.32 GeV/c 2 is seen. The shaded histogram represents the 
events in the — > K + K~it + mass sidebands combined with the ir° candidates. In 
Fig. the mass distributions result from the combination of the candidates 
with the photon pairs from the tt° signal and sideband regions of Fig. IHb) (the side- 
band distribution is again shaded). In Figs, QJc) and^d) the 2.32 GeV/c 2 signal 
is absent from the sideband distributions indicating quite clearly that the peak is 
associated with the D+ir system. In order to improve mass resolution, the nominal 
.D+ mass [2] has been used to calculate the energy. 

The -D+7T mass distribution for p*(D+7r°) > 3.5 GeV/cis shown in Fig.|2fa). 
The fit function drawn on Fig. |2f a) comprises a Gaussian function describing the 
2.32 GeV/c 2 signal and a polynomial background distribution function. The fit 
1 The inclusion of charge-conjugate configurations is implied throughout this paper. 
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Figure 1: BABAR experiment, (a) The distribution of K + K~ir + mass for all can- 
didate events. Additional selection criteria have been used to produce the lower 
histogram, (b) The two-photon mass distribution from D+7r° candidate events. Df 
and 7T° signal and sideband regions are shaded, (c) The Dfir mass distribution for 
candidates in the Df signal (top histogram) and K + K~n + sideband regions (shaded 
histogram) of (a), (d) The -D^77 mass distribution for signal candidates and 
a photon pair from the 7r° signal region of (b) (top histogram) and the sideband 
regions of (b) (shaded histogram). 



yields 1267 ±53 candidates in the signal Gaussian with mass (2316.8 ±0.4) MeV/c 2 
and standard deviation (8.6 ± 0.4) MeV/c 2 (statistical errors only). The signal, la- 
belled as .D*j(2317) + , is observed in both the 07r + and K*°K + decay modes of the 
Df. In addition, a sample of — > K + K^7i + tt° decays is selected by adding n° 
candidates to each K + K~ir + candidate. Each resulting candidate is combined 
with a second 7r° candidate with lab momentum greater than 300 MeV/c. A clear 
Z^*j(2317) + signal is observed as shown in Fig. Efb). A Gaussian fit yields 273 ± 33 
events with a mean of (2317.6 ± 1.3) MeV/c 2 and width (8.8 ± 1.1) MeV/c 2 (statis- 
tical errors only). The mean and width are consistent with the values obtained for 
the .D+ — > K + K~tt + decay mode. The mass distribution of the — > K + K~7t + tt° 
sample (not shown) peaks at (1967.4 ± 0.2) MeV/c 2 (statistical error only). 
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Figure 2: BAR4R experiment. The Dfn° mass distribution for (a) the decay D s 
K+K-1T+ and (b) the decay D+ -> K+R-n+ir . 



Monte Carlo simulations have been used to investigate the possibility that 
the -D*j(2317) + signal could be due to reflection from other charmed states. This 
simulation includes e + e~ — > cc events and all known charm states and decays. The 
generated events were processed by a detailed detector simulation and subjected to 
the same reconstruction and event-selection procedure as that used for the data. No 
peak is found in the 2.32 GeV/c 2 Dfn signal region. 

Mass resolution estimates for the K + K~-k + -k system are obtained directly 
from the data using a fit to the mass distribution — > K + K~7t + tt°. The measured 
width from this mode is consistent with that of the _D*j(2317) + signal. It can be 
concluded that the intrinsic width of the D* J (2317) + is small (T < 10 MeV). 

A search has also been performed for the decay _D*j(2317) + — > Dfj. 
Shown in Fig. Efa) is the -D^7 mass distribution obtained by combining a candi- 
date in the signal region of Fig. E^a) with a photon with an energy of at least 150 MeV 
that does not belong to a 77 combination in the signal region of Fig. ^b). The re- 
quirement that the p* of the -D+7 system be greater than 3.5 GeV/c is also imposed. 
There is a clear D*(2112) + signal, but no indication of -D*j(2317) + production. The 
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Figure 3: BABAR experiment. The mass distribution for (a) Df'y and (b) -Dj~77 
after excluding photons from the signal region of Fig. IHb). (c) The D^tt ^ mass 
distribution. The lower histograms of (b) and (c) correspond to -D+ 7 masses that fall 
in the D*(2112) + signal region as described in the text. The vertical line indicates 
the ££,(2317)+ mass. 

-D+77 mass distribution for p*(D+77) > 3.5 GeV/c, excluding any photon that be- 
longs to the 7T° signal region of Fig. db), is shown as the upper histogram of Fig.Efb). 
No signal is observed near 2.32 GeV/c 2 . The shaded histogram corresponds to the 
subset of combinations for which either -D^7 combination lies in the D*(2112) + re- 
gion, defined as 2.096 < < 2.128 GeV/c 2 . Again, no ^(2317)+ signal is 
evident, thus demonstrating the absence of a Z)*(2112) + 7 decay mode at the present 
level of statistics. The D+n -/ mass distribution, excluding any photon that belongs 
to any 7r° candidate, is shown as the upper histogram of Fig. Efc). The shaded 
histogram corresponds to the subset of combinations in which the -D^7 mass falls 
in the D*(2112) + region. No signal is observed near 2.32 GeV/c 2 in either case. A 
small peak, however, is visible near a mass of 2.46 GeV/c 2 . This mass corresponds to 
the overlap region of the D*(2112)+ -> £+7 and £^(2317) + -> D+n signal bands 
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that, because of the small widths of both the L>*(2112) + and _D*j(2317) + mesons, 
produces a narrow peak in the D+7r 0/ y mass distribution that survives a D*(2112) + 
selection. 

If the peak in the Dfn ^ mass distribution of Fig. [3(c) were due to the pro- 
duction of a narrow state with mass near 2.46 GeV/c 2 decaying to D*(2112) + 7r°, the 
kinematics are such that a peak would be produced in the D+n mass distribution 
at a mass near 2.32 GeV/c 2 . Such a D+n mass peak, however, would have a root- 
mean-square of ~ 15 MeV/c 2 , which is significantly larger than that obtained for the 
_D*j(2317) + signal. In addition, Monte Carlo studies indicate that if the apparent 
signal at 2.46 GeV/c 2 were due to a state that decays entirely to D*(2112) + 7r°, it 
would produce only one-sixth of the observed signal at 2.32 GeV/c 2 . 

The BABAR experiment was somewhat coutious in claiming the discovery 
of this second state. "Although we rule out the decay of a state of mass 2.46 GeV/c 2 
as the sole source of the D+7r° mass peak corresponding to the -D*j(2317) + , such 
a state may be produced in addition to the -D*j(2317) + . However, the complexity 
of the overlapping kinematics of the D s *(2112)+ -> D+7 and D* J (2317)+ -> D+tt° 
decays requires more detailed study, currently underway, in order to arrive at a 
definitive conclusion pQ." 

4 Confirmation of D* sJ (2317) + and observation of D sJ (2A58) + by other 
experiments. 

Using 13.5 and 78 fb _1 integrated luminosity, the CLEO and Belle ^1] experi- 
ments readily confirmed the existence of D*j(2317) + (shown in figEJ). I n terms of 
Am = m(K + K-n + n°)-m(K + K-n + ), CLEO reported a value of Am = 350.3 ±1.0 
MeV/c 2 with 231 ± 30 events. Belle reports Am = 348.9 ± 0.5 MeV/c 2 with 643 
± 50 events. These values are in good agreement with the BABAR measurement of 
Am = 384.4 ± 0.4 MeV/c 2 . 

In addition, both the CLEO and Belle Collaborations have analyzed the 
D*(2112) + tt° mass distribution finding evidence for structure in the 2.46 GeV region 
(see fig. EJ). Defining now Am = m(K + K~7r + 7r 0r y) — m(K + K~7r +r y), they reported 
the following parameters for this state: Am=351.2 ± 1.7 MeV/c 2 with 41±12 events 
(CLEO) and Am=344.1 ± 1.3 MeV/c 2 with 79 ± 18 events (Belle). 

Belle also reported evidence for both states ^Uj m the B decays to B — > 
DD*j(2317) + and B — > DD s j(2A58) + . In addition, they reported an observation 
of .D s j(2458) + — > -D^7 in both B decays and continuum e + e~ annihilations. The 
angular analysis shows the expected behaviour for a J p = 1 + state. 
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Figure 5: Am(D* + 7r°) from CLEO and Belle Collaborations showing the evidence 
for D s j(2458) + . Data are from continuum e + e~ annihilations. The distribution 
from D* + sidebands is also shown. 
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No evidence has been found for narrow structures in or D+irir final 

states. 

5 The observation of _D s j (2458) + by the BABAR experiment. 

In order to study the -D+7r°7 system, each candidate is combined with all com- 
binations of accompanying tt° candidates with momentum greater than 300 MeV/ c 
and photon candidates of energy greater than 100 MeV/c 2 . To suppress background, 
photon candidates that belong to any n° candidate are excluded and it is required 
that the combined momentum p* in the e + e~ center-of-mass mass frame of each 
D+7r°7 combination be greater than 3.5 GeV/c. The -D+7r°7 invariant mass dis- 
tribution is shown in Fig. El A small peak is observed near 2.46 GeV/c 2 . The 
background underneath this peak is from several sources, which can be described in 
terms of mass differences defined as: 

Am(L>+7) = m(K-K + 7r + j) - m(K~K + -K + ) (1) 
Am(D* s + TT°) = m(K-K + Tr + ir ^) - m(K~ K + ir + ^) . (2) 

A scatter plot of these quantities is plotted in Fig.^b. Two signals are clearly visible: 
.D*(2112) + — ► -0^7 decay combined with unassociated 7r° candidates, which appears 
as a horizontal band, and _D*j(2317) + — > D^tt decay combined with unassociated 
7 candidates, which appears as a band that is almost vertical. 

An enhancement is evident in the vicinity of the overlap of these two bands 
with a D+7r°7 mass near 2.46 GeV/c 2 . The upper histogram of Fig. Eh shows events 
in the D*(2112) + signal region, and the shaded histogram shows those in the two 
D*(2112) + sidebands. One can conclude that a state with decay to D+7r°7 is seen 
to exist over a background from D*j(2317) + combined with a 7. This background 
is peaked at a mass slightly higher than the signal. In Fig. Ell the subtracted plot 
shows the narrow signal fitted to a Gaussian shape on a second order polynomial 
background at Am(D+7) = 346.2 ± 0.9 MeV/c 2 (statistical errors only). 

The signal, which is labelled _D s j(2458) + , may decay either through _D+7r°7, 
D*(2112) + 7T° or _D*j(2317) + 7. To disentangle these decay modes and extract the 
most information from the data, a binless-maximum likelihood fit to the D+7r°7 
system has been developed. The fit determines a _D s j(2458) + mass of (2458.0 ± 
1.0) MeV/c 2 and a measured Gaussian width equal to (8.5 ±1.0) MeV/c 2 . This mass 
value (yielding Am = 345.6 ± 1.0 MeV/c 2 ) agrees with that obtained by BELLE 
but differs from that obtained by CLEO. 
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Figure 6: BABAR experiment, (a) The mass distribution for all selected .D+7r°7 
combinations. The shaded region is from sidebands defined by 1.912 < 
m(K-K + n + ) < 1.933, 1.999 < m(K~K + n + ) < 2.020 GeV/c 2 . (b) The value of 
Am(Df'-f) versus Am(D* + n°) for all combinations. The horizontal lines delineate 
three ranges in Am(Df'-f). (c) The Am(D* + 7r°) mass distribution for the middle 
range of Am(f s + 7) (white) and for the upper and lower ranges (shaded), (d) The 
difference of the two histograms shown in (c). The curve is the fit described in the 
text. 



The shape of the D+ir and -D^7 distributions in the -D s j(2458) + mass re- 
gion can be used in order to distinguish the two possible D s j(2A58) + — > D*(2112) + tt° 
and _D s j(2458) + — * -D*j(2317) + 7 decay modes. These shapes are influenced by the 
allowed kinematic range for D s j(2458) + decay, as shown in Fig.[7^. Figs. [7|D-[7t: show 
the _D s j(2458) + background subtracted D^tt and Dfj mass projections compared 
with MC simulations of the two hypotheses. The decay D sJ (2A58) + -> £>*(2112) + 7r° 
is clearly favored. 

The distribution of the angle $h of the decay _D*(2112) + — > in its 

center-of-mass with respect to the -D s j (2458) + can be used investigate the spin- 
parity of the _D s j(2458) + . The resulting efficiency-corrected cos-dh distribution is 
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Figure 7: BABAR experiment, (a) The D+n versus -D+7 mass distribution 
for all D^7r°7 combinations. The curves indicate the kinematically allowed re- 
gion for _D s j(2458) + decay. (b) Sideband subtracted Dfj mass distribution 
with (line) Monte Carlo simulation for D sJ (2458)+ -> £>*(2112) + 7r° and (dashed) 
£> s j(2458) + -> D;j(2317) + 7. (c) A similar plot for the £>+7r° mass distribution, (d) 
The efficiency corrected yield as a function of dh (statistical errors only). The solid 
(dashed) histogram corresponds to the best fit of a sin 2 -& h (1 + cos 2 $/J distribution. 



shown in Fig. 01 (statistical errors only). This distribution is not consistent with a 
sin^ distribution, which rules out a -D s j(2458) + spin-parity assignment of J p = 0~ . 

MC simulations have been used to measure the detector resolution, lead- 
ing to the conclusion that the intrinsic width of the D s j (2458) + is small (r < 
10 MeV/c 2 ). 

6 Conclusions 

The decay of any cs state to Dfir or D*(2112) + 7r° violates isospin conservation, 
thus guaranteeing a small width for these states. It is possible that the decays 
proceeds via t] — tt° mixing, as discussed by Cho and Wise ■ 
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The low mass and the absence of a D* sJ (2317) + -> £+7 favors J p = 0+ 
for D*j(2317) + . The mass of D*j(2317) + lies below the DK threshold, the mass of 
D sJ (2A58) + lies above DK and below D*K thresholds. 

Different interpretations for these states have been proposed. In ref. 113*1 114*] 
models in terms of baryonia or molecules have been proposed. Ref. [!""] provides an 
explanation in terms of relativistic vector and scalar exchange forces. Ref. [H"] uses 
HQET plus chiral symmetry to predict parity doubling, i.e. the expectation is that 
the mass splitting between the 1 + and 1~ states should be the same as for the + 
and 0~ states. 

Since a cs mesons with these masses contradicts current models of charm 
meson spectroscopy El HI , most likely these models need modification. 
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